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ABSTRACT 

We compare the surface brightness-inclination relation for a sample of COSMOS pure disk galaxies 
at 2^0.7 with an artificially redshifted sample of SDSS disks well matched to the COSMOS sample 
in terms of rest-frame photometry and morphology, as well as their selection and analysis. 
The offset between the average surface brightness of face-on and edge-on disks in the redshifted SDSS 
sample matches that predicted by measurements of the optical depth of galactic disks in the nearby 
universe. In contrast, large disks at z ~ 0.7 have a virtually flat surface brightness-inclination relation, 
suggesting that they are more opaque than their local counterparts. This could be explained by 
either an increased amount of optically thick material in disks at higher redshift, or a different spatial 
distribution of the dust. 

Subject headings: cosmology: observations — galaxies: spiral — galaxies: evolution — galaxies: ISM 
— dust, extinction — surveys 



1. INTRODUCTION 

After a long-lasting controversy o n whether galactic 
disks are optically th ick or thin (e.g. iDisnev et al.lll989l : 
iBurstein etall 119911) . a consensus emerged during the 
last decade that disks in the local universe (i) behave 
like optically thick systems as far as integrated pho- 
tometric properties in the ultraviolet (UV) and visual 
are concerned (e.g . iShao et al] 120071: iDriver et all 120071: 
iMaller et al.ll2009h . while at the same time (ii) periph- 
eral and inter-arm regions in spiral galaxies are more 
transparent than spiral arms or the core of the disk (e.g. 
iWhite et all 120001 : IHolwerda et al.1 120051 and references 
therein) . 

Beyond the nearby universe, observational constraints on 
the extinction in disk gal axies haye bee n hard to obtain 
(see, e.g., the review bv lCalzettil[200l as (i) some of 
the locally employed measurement techniques can only 
be gainfully applied over a limited redshift range (e.g. 
IHolwerda et al.l [20071) : (ii) the spectral energy distri- 
bution (SED) of distant galaxies between the UV and 
infrared (IR) is often only partially sampled; and (Hi) 
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statistical tests - such as attenuation-inclination rela- 
tions - could onl y be applied to fairly small samples (e.g. 
ILillv et al.lfl998l) . 

Pan-chromatic wide field surveys like COSMOS in princi- 
ple provide both the necessary multi- wavelength informa- 
tion, as well as the means to select a statistically signifi- 
cant and (morphologically) well-defined sample of disk 
galaxies thanks to high-resolution HST (Hubble Space 
Telescope) imaging. Here we perform a comparative 
study of the surface brightness-inclination relation in a 
sample of local and distant (z~0.7) pure disk galax- 
ies with the aim of testing if their average opacity has 
evolved in the past 6.5 Gyr. Our galaxy samples are de- 
scribed in §[3 followed by results (§[3]) and conclusions 

We adopt a concordance cosmology with fi m =0.25, 
fiA+fim = 1 and Hq = 70 k m s~ 1 Mpc~ 1 . Magnitudes are 
given in the AB-system of lOkd (j 19 741) . 

2. THE DISK SAMPLES 
2.1. Strategy 

A study of the opacity of distant disk galaxies that re- 
lies on structural information (e.g., surface brightness) 
requires both high-resolution imaging and a consistent 
sampling of rest-frame wavelengths. Avoiding so-called 
'morphological X-corrections' is especially imp ortant in 
view of, e.g., size gradients with color (e.g., Ide Jond 
|1996| ). and also because unrecognized bulge components 
can ind uce spurious signa ls in attenuation-inclination re- 
lations ([Tuffs et al.H200l . 

The limitation of COSM OS HST imaging to a single 
band (the F814W filter; see lScoville et al.ll2007t) prevents 
the construction of disk galaxy samples over a contin- 
uous redshift range. We hence resort to a comparative 
study, taking advantage of the fact that the central wave- 
length of the SDSS g-band matches the rest-frame wave- 
length of objects observed in the F814W filter at red- 
shift z^0. 7. By artificially redsh ifting local galaxies to 
z ~ 0.7 (cf. iKampczvk et al.|[2007l ) it becomes possible to 
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directly compare the attenuation properties of distant 
disks to those in the nearby universe. Moreover, since 
our COSMOS disks and the redshifted local sample are 
well matched in terms of both photometric and morpho- 
logical -ftT-corrections, comparisons between the two data 
sets should eliminate systematic biases. 

2.2. COSMOS Disks at 0.6<z<0.8 
2.2.1. Morphological Measurements 

The z~ 0.7 COSMOS disks used for the present analy- 
sis are selected from a complete sample of 55,324 galaxies 
with / < 22.5 listed in the Zurich Structure and Mor- 
phology Catalog vl.O ('ZHSM catalog'; Carollo et al., 
in prep.). Morphological measurements are carried out 
on HST/ AGS F814W (/-band) images with a pixel scale 
of 0. 05" /px and resolution of ~0.1" (jKoekemoer et all 
120071 ). 99.5% of all catalog entries have a classification 
as either 'early-type', 'late-type' or 'irregular/peculiar' 
according to the ZEST algorithm (Zurich E stimator of 
Structural Types (vl); iScarlata et all 120071) . Briefly, 
ZEST performs a principle component (PC) analysis on 
five non-parametric structural estimators: asymmetry A, 
concentration C, Gini coefficient G, second-order mo- 
ment of the brightest 20% of the pixels associated with 
an individual galaxy, M20, and cllipticity e. Late- type 
galaxies (ZEST CLASS = 2) arc further divided into four 
bulgcncss categories based on the distribution of Sersic 
indices in the different regions of a three-dimensional PC 
spac43 occupied by these galaxies. 

The compl ete sample ha s also been modeled with single- 
com ponent iSersid (1 1968T) profi l es using the Gim2d pack- 
age (|MMliauXlimarj[l99l [UmardiFin EOQl) . We 
thus obtain total flux i*totj semi-major axis half-light ra- 
dius -R1/2; ellipticity e = l — b/a (where a and b are the 
semi-major and semi-minor axes of the fitted Sersic pro- 
file); and the value of the Sersic index n (for details see 
ISargent et al.l[20 07) . For the current analysis we will only 
use the recovered size and axis ratio, which are essential 
to estimating the surface brightness and orientation of 
our disks. 

We select all COSMOS disk galaxies with a negligi- 
bly small bulge component (ZEST CLASS = 2.3; median 
Sersic index n< 1) that are not candidates for being ei- 
ther spurious sources (flags JUNK > and ACS_CLEAN = 0) 
or stars (flag STELLARITY> 0). This results in a sample 
of 6015 pure disks with / < 22.5 that are not located in an 
area where the COSMOS multi-wavelength photometry 
(cf. § I2.2.2[) is affected by charge bleeds from saturated 
stars. 

2.2.2. Multi-Wavelength Photometry & Redshift Estimates 

The i + -band-sel ected COSMOS photometry catalog 
(jllbert et al.l 120091 : P. Capak et al. 2010, in prep.) 
contains >500,000 sources with i < 25, is complete 
down to this magnitude and provides PSF-matched 
(FWHM = 1.5") photometry in 30 broad, medium and 
narrow band filters. The wavelength range covered 
by these observations extends from 1550 A to 8/xm 

12 The ZEST classification scheme reduces the initially five- 
dimensional space of structural parameters to three dimensions. 
The axes defining the new metric are linear combinations of A, 
C, G, M20 and e, and correspond to the three dominant principle 
components. They represent 92% of the variance in the data set. 



(iTaniguchi et al.ll2007t iCapak et aT]|2007t iSanders et al.l 
2007). 

We cross-correlated the ZHSM catalog with the photom- 
etry catalog in order to obtain multi- wavelength and dis- 
tance information for our disk galaxy sample. Straight- 
forward positional matching using a search radius of 
0.6" resulted in ground-based optical identifications for 
1,598 pure disk galaxies in the redshift range of inter- 
est, 0.6 <z< 0.8. Twenty-three of these were excluded 
because the association between ACS and ground-based 
candidates was not unique. Most of these ambiguous 
matches consist of two ACS objects being blended into a 
single source in the ground-based imaging. 
Successfully matched disks galaxies generally have a 
photometric r edshi ft estimated with the code Le Phare 
(jllbert et all [2009). The accuracy of the photo- 
metric redshifts was calibrated with more than 4000 
high-co nfidence spectrosc opic redshifts from zCOSMOS 
sources (jLillv et al.ll2009D with i + < 22 .5. In this magni- 
tude range and at z < 1.25. lllbert et al.l found that photo- 
z measurements have a dispersion of a( Az /{i+z)) ~ 0.007. 
Based on this statistically expected accuracy of the pho- 
tometric redshifts, we exclude objects with uncertain 
photo-zs. To remain in the sample, the width of the 
photo- z probability distribution for a specific source must 
be <2(T of the dispersion. 30 (~2%) of the 1575 disks 
were rejected due to this criterion, leaving 1545 pure 
disk galaxies with good quality photometric redshifts 
0.6<2<0.8. 

2.3. SDSS Reference Sample 

iKampczvk et al.l (|2007t ) use SDSS g-band images of 
a volume-limited sample of local disk galaxies (median 
redshift z ~ 0.02) to simulate the appearance of these 
galaxies in COSMOS ACS /-band images if redshifted 
to z-0.7. As mentioned in §CD]the ACS F814W filter 
almost perfectly samples th e rest-frame SDSS q-ban d at 
this redshift. Furthermore, IKampczvk et alj ( 20071 ) ac- 
count for effects of band-pass shifting, cosmological sur- 
face brightness dimming, changing spatial resolution and 
increased noise. The simulations used here do not in- 
clude luminosity evolution as its implementation is not 
straightforward (e.g. non-uniform brightening of galaxy 
pixels as conceivably caused by a dust component of un- 
known geometry and composition). 
The SDSS z= o.7 galaxies were processed in exactly the 
same way as the real COSMOS data. This not only 
applies to structural measurements and fits to the 
SEP, but also to sour ce detection with SExtractor 
(|Bertin fc Arnouti |1996[ ) following the insertion of the 
simulated galaxies into COSMOS ACS pointings. In 
particular, the same SExtractor configuration parame- 
ters used to generate the original list of COSMOS ACS 
detections in iLeauthaud et all (|2007l ) were adopted for 
the extraction of the SDSS z= o.7 objects. About 6% of 
the SDSS z= o.7 galaxies were not recovered by SExtrac- 
tor. The success rate of the morphological measurements 
was 100%. We thus obtain a comparison sample of 75 
pure disk galaxies that, after redshifting, are brighter 
than 7 = 22.5. 

3. RESULTS 

In the following we consider the surface brightness- 
inclination relation in the rest-frame _B-band. This has 
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the twofold advantage of bordering on the SED region 
(g-band) sampled by our SDSS and COSMOS imaging, 
and of providing results that are readily comparable with 
predictions and previous findings in the literature. 
The average surface brightness is obtained by distribut- 
ing half (hence the second term in eq. H]) of the total 
/-band flux over an ellipse with a semi-major axis R\n 
(in arcseconds). The semi- minor axis follows from the 
Gim2d value of the axis ratio: 



fx(B) = I + 2.5 log (2) + 2.5 log ^Rf /2 ) (1) 

+ 2.5 log (b/a) - 10 log (1 + z) - K F8UW! b (z) . 

The /C-correction Kf8hw,b(z) is determined using the 
best -fit SED in the templat e library of the ZEBRA pack- 
age l)Feldmann et al.l 12006ft . The relation between p(B) 
and that of the system when viewed face-on (//°) is often 
parametrized as: 



H(B) = n(cos(ij\ = fi fo + C x 2.51og(cos(i) 



(2) 



Here i is the inclination angle (i = for face-on and 
*/2 for edge-on systems). The amount of attenuation 
is parametrized by the opacity constant C which ranges 
from unity for transparent disks to zero for opaque disks. 
For a non-zero C, the apparent surface brightness of an 
inclined disk is thus always more intense than in the face- 
on object. 

The distribution of axis r atios b/a can be a n indication of 
the opacity of disks fe.g. Uones et al.lfl996f) . In Fig. [Tjwc 
plot the axis ratios for our COSMOS disks (right) and 
the SDSS z= o.7 sample (left). Instances o f b/q~ l and 
are rare due to intrinsic disk ellipticity (jRvdenl 120041 ) 
and the finite thickness of edge-on systems. A one-sided 
Kolmogorov-Smirnov (K-S) test indicates that - within 
the errors shown - there is a non-negligible probabil- 
ity that the measured axis ratios in the entire samples 
are drawn from a flat distribution (SDSS z= o.7: p = 0.19; 
COSMOS: p = 0.16), as illustrated by the cumulative dis- 
tribution function crossing both panels on an approxi- 
mate diagonal. Fig. Q] also shows the histograms (in 
grey) of b/a for the subset of disks with physical half- 
light radii r 1 / 2 >5kpc. We introduce this restriction in 
order to obtain a complet^^ sample of disks for the fi- 
nal analysis. While the axis ratio distribution of large 
SDSS z= o.7 disks remains flat (K-S probability p = 0.99), 
that of the corresponding COSMOS population is skewed 
toward elongated objects (p < 10~ 5 for the distribution 
being flat). The different behaviour hints at opacity vari- 
ations which we quantify in the next paragraph. 
Axis ratios are related to inc lination ( cf. als o scale along 
upper edge of Fig. [Q) by the lHubbld (|1926ft formula: 



'(0 



(b/a)' - (b/a) 



2 

min 



1 - (b/a) 



2 

min 



(3) 



with (b/a) m i n = 0.15, in keeping with measure ments of 
the minimal flattening in late-type spirals (e.g. I Guthrie! 

13 The threshold of 5kpc corresponds to the minimal size 
at z~0.7 which ensures that the selection limit in magnitude 
(7 = 22.5) can be reached for the entire range of 'normal' surface 
brightness valu es (i.e. negl e cting likel y rare very low s urface bright- 
ness disks), cf. ILillv et all 1 (1993 ) and lSargent et all l(2007h . 



11993 I^IaFFzEai200l . 

We present the surface brightness- inclination relation for 
large SDSS z= o. 7 and COSMOS disks (51 and 611 objects, 
respectively) in Fig. [2l The error bars on the median 
surface brightness values in bins of increasing inclination 
(cf. also Table QJ span the 95% confidence region, esti- 
mated with N cos ^ x 100 bootstrap realizations (N cos ^ 
is the number of disks in a given bin of inclination). 
The uncertainties on the COSMOS median are signifi- 
cantly smaller, in agreement with the expectation that 
they should scale roughly as 1 /y / N cos(i) . 
The dashed blue line in Fig. [5] (left) marks the sur- 
face brightness-inclination relation for a B-band cen- 
tral face-on optical depth of = 3.8, the average opac- 
ity of nearby disks recently determined by iDriver et al. | 
(12007ft . The plotted l ine is taken from the predictions 1 ^1 
in IMollenhoff et"aT1 (|2006D . based on radiative trans- 
fer models of st able local disks (jPopescu et al.1 120001 : 
iTuffs et al.1 12004ft . The good agreement with our mea- 
surements demonstrates that we correctly recover the 
attenuation properties of local disks even when they are 
redshifted to z^0.7. 

We now determine the average correction between face- 
on and observed surface brightness. At first sight, Fig. 
[5] might suggest that - within uncertainties - the surface 
brightness-inclination relation is flat both for SDSS z= o.7 
and COSMOS disks (see dotted horizontal line marking 
the global surface brightness average). However, as we 
will now show, the opacity constants describing local and 
z~0.7 disks turn out to be significantly different. 
To properly account for strongly asymmetric error bars 
(see Fig. [5]) we interpret the distribution of bootstrapped 
medians as a probability distribution and accordingly 
draw pairs of medians (cos(i), ~p(B)) in all bins of in- 
clination. We then fit eq. © to each of the sets of 
resampled medians and subsequently locate the peak of 
the resulting distributions of the free parameters C and 
/i fo . 

Fig. [3] (main window) shows the best-fitting parameter 
pairs ((n io (B)} 7 C) from 10,000 realizations. The cores of 
the point clouds are highlighted with (smoothed) equal 
density contours that enclose 68, 90 and 95% of the 
points. By projecting the scatter plots along the ver- 
tical and horizontal axes we obtain 95% confidence in- 
terval^ for the free parameters (//°(i?)) and C for the 
local and high-z disk samples (Table [3]). The average 
surface brightness of face-on disks, (p, (B)) , increases 
by ~ 1 mag between z ~ and 0.7 (cf. also ILillv et all 
119981 : iBarden et al.l 120051 ). For the opacity constant we 
find C(z ~ 0.7) G [0.18,0.62] and C(z ~ 0) G [0.01, 0.13] 
and for the most probable valuet03 C(z~0) = 0.47 and 

14 The fact that our pure disks are fit by a single Sersic com- 
ponent implies that average and central surface brightness (j.q(B ) 
are related by a con stant offset (see, e . g.JGraham fc Driverll2005D . 
The predictions of IMollenhoff ct~aT1 1 12009) for the inclination- 
dependence of ~po(B) can thus be meaningfully compared to our 
findings. 

5 Note that the 95% confidence limits include the diffuse cloud 
of outliers at (fi to (B)) « 21.8 (Fig. 0, occurring as a result of small 
number statistics in the first and last SDSS z= o.7 inclination bin. 

16 The evolution of the opacity constant C is not brought about 
by differing distributions of axis ratios in the two samples. Resam- 
pling the surface brightness values measured for large COSMOS 
disks according to the distribution of b/a in the SDSS z= o.7 sam- 
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C(z ~ 0.7) = 0.07. The corresponding surface brightness- 
inclination dependence (eq. [2]) is indicated in red in Fig. 

m 

We note that by usin g the 'unbrightcned' sample of 
IKampczvk et al.l (|2007D (see § !2.3p only the most lumi- 
nous local disks satisfy the selection criterion 7<22.5 
when redshifted to z~ 0.7. The general l uminosity evo- 
lution of the disk galaxy population (e.g. iScarlata et al.l 
120071) implies that only the brightest COSMOS disks 
at 0.6 < z < 0.8 have masses similar to those of the 
SDSS z= o.7 sample. As the degree of extinction de- 
pends on the am ount of attenuating material (e.g. 
iMasters et al.ll20Toh . this could be partly responsible for 
the derived difference in the opacity constant C. We can 
generate a sample of COSMOS disks that, to first order, 
have masses similar to the SDSS 2= o.7 galaxies by bright- 
ening the selection limit for COSMOS galaxies by ^Az 
magnitudes to 7w21.8. This restriction does not change 
our finding of a nearly fiat (C(z - 0.7) = 0.05±g;?|, in this 
case) surface brightness-inclination relation at z~ 0.7. 



4. CONCLUSIONS 

We have investigated the average opacity of distant 
(0.6 < z < 0.8) COSMOS spiral galaxies by direct com- 
parison with local disks artificially redshifted to z~0.7. 
The outcome of inclination-based attenuat ion measure- 
ment s is susceptib l e to s election effects (e.g. lDavies et al.l 
Il993t Uones et"al] 119961) . By processing our two sam- 
ples identically and by applying a morphologically and 
photometrically consistent sample selection, we have 
taken the necessary precautions to reduce such system- 
atic biases. We recover the expected surface brightness- 
inclination relation for nearby disks even after the simu- 
lated shift to z^0.7 (in the parameterization of cq. © 
this implies an opacity constant C(z ~ 0) = 0.47+Q2g). 
For high-z COSMOS disks a significantly lower value of 
C(z~ 0.7) = 0.07±0.06 is found, implying, on average, a 
nearly constant relation between surface brightness and 
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inclination as is expected for optically thick spiral galax- 
ies. 

Previous studies suggest that the extinction laws in star- 
forming galaxies at similar redshifts as our COSMOS 
sample do not differ str on gly from thos e in local sys- 
tems (e.g.. iCon rov 2009 or Calzetti 200 1] and references 
therein) . It thus seems unlikely that the increased opac- 
ity of our z-0.7 COSMOS disks is due to a different 
chemical composition of the dust. Given that our low- 
and high-z disk samples have similar luminosities we can 
also rule out that the evolution is due to the locally ob- 
served scali ng of dust opacity with the blue luminosity 
of galaxies (|Wang k, Heckmanlll996l) . 
Other possible explanations for the flat surface 
brightness-inclination relation at z ~ 0.7 are the presence 
of more attenuating material, or a different spatial ar- 
rangement thereof. Evidence t h at bo th factors might 
contribute exists: iGenzel et al.1 ()2008l ) report stronger 
turbulent motion in disk-like systems at high redshift 
that could increase the scale height of the dust, and re- 
cent measurements of molecular line emission in typical 
late type galaxies at z ~ 1.5 have revealed l arge gas frac- 
tions in excess of 50% of the baryonic mass (jDaddi et al.l 
120091 ). The measurements presented here are not suffi- 
cient to infer the relative importance of such potential 
contributions. Additional constraints from complemen- 
tary measurements or different wavelength regions are 
thus indispensable to determine the causes for the opac- 
ity evolution we observe between z~ and 0.7. 
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TABLE 1 

Median values of inclination (col. 2) and 
OF THE ASSOCIATED average rest-frame 
B-band surface brightness (col. 3; CF. 
also Fig.[2J. Errors span the 95% 
confidence interval. 



samplc/rcdshift 1 — cos(i) ^{B) 

SDSS Z=0 .7 0.148±°°«2 22.2371°,-°?? 

0.310™ 21.866lg;?l? 

n ctri +0.034 91 79 „+0.213 

u.ooi_ 037 zi.(zo_ 307 

793+O.O58 (j-j ^qtj+0.393 



COSMOS 

(*~0.7) 4381°°1b 

,+0.013 
3 -0.018 

0.746+nnna 21.042 



, +0.138 



-0.019 

38l , 
0.603l 0018 

+0.007 r,-, n/i r>+0.082 

0.009 z±iU * z -0.070 

0.8831°;°}* 21.1091°;™ 



TABLE 2 

Favoured values (cf. Fig. [3j of the free 
parameters {fj° (b)} and c in the surface 

brightness-inclination relation of eq. j2j. 

Errors span the 95% confidence interval. 



sample /rcdshift 


(M fo (B)) 


C 


SDSS z =o.7 


oo i 7 o+0.158 
zz.kz_ 339 


465+ ' 158 

U.4DO_q 2gl 


COSMOS (z~0.7) 


91 1fi«+ 091 

2l.lb8_ 0g6 


n n65+ 060 

U.UDO_ Q ogl 



1 - cos(i) 




b/a 

Fig. 1. — Distribution of axis ratios b/a for pure disk galaxies at 0.6<z<0.8 in the COSMOS field {right) and SDSS z= o.7 disks {left). 
The histograms are normalized to the total number of galaxies in the two samples. The distributions of axis ratios of disks with ri t% > 5 kpc 
are indicated in grey. The upper axis shows the correspondence between inclination i and axis ratio b/a. 

The cumulative distribution function (CDF) of axis ratios in the total disk population is plotted with a solid black line (see scale at right). 
Dashed lines delimit the associated 99% confidence region. 
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b/a 



to 




1 — cos(i) 

Fig. 2. — Rest frame B-band surface brightness-inclination relation for large (ri/2 ^ 5kpc) artificially redshifted local disk galaxies (left) 
and observed large disks at z~0.7 (right). Dotted lines: median surface brightness averaged over all inclinations. Black circles: median 
surface brightness and inclination in discrete bins of inclination (error bars span the 95% confidence interval). Red lines represent the 
best-fitting surface brightness-inclination relation in eq. ||2) for local and z ^ 0.7 disks. On the left, the blue line shows the prediction 
by iMol lcnhoff ct al. ( 2006) based on the mean B-band opacity of local spiral disks in the Millenium Galaxy Catalogue (t|^ = 3.8±0.7; 
IDriver et a!ll2l)07T) . 



a 

CO 

M 

a 

o 

o 



o 

CO 

o 



0.1 

0.05 


1 



0.5 - 



— i — m— 



T-t- i I I 




• 0.6 < z < 0.8 (COSMOS) 




I 1 1 1 1 I 1 1 1 1 I 



22.5 



22 



21.5 




<M f °(B)> 



0.05 0.1 



fraction 



Fig. 3. — Distribution of the best-fit parameter pairs ({/x fo (B)), C) in eq. J2t . obtained by repeated fits to the resamplcd data in Table[T] 
(see text for details). Light and dark grey is used for the low and high redshift sample, respectively, as in Fig. [2] The contours in the main 
window arc isopycnics enclosing 95, 90 and 68% of the points. The panels along the edges of the figure show the projected distributions of 
average face-on surface brightness in the B-band, (fi fo (B)) (top), and of the opacity constant C (right). 



